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We study smectic liquid crystal order in a cell with a heterogeneous substrate imposing surface 
random positional and orientational pinnings. Proposing a minimal random elastic model, we 
demonstrate that for a thick cell the smectic state without a rubbed substrate is always unstable 
at long scales, and for weak random pinning is replaced by a smectic glass state. We compute 
the statistics of the associated substrate-driven distortions and the characteristic smectic domain 
size on the heterogeneous substrate and in the bulk. We find that for weak disorder, the system 
exhibits a three-dimensional temperature controlled phase transition between a weakly and strongly 
pinned smectic glass states akin to the Cardy-Ostlund phase transition. We explore experimental 
implications of the predicted phenomenology and suggest that it provides a plausible explanation 
for the experimental observations on polarized light microscopy and x-ray scattering. 

PACS numbers: 64.70.pp, 61.30.Hn, 64.60.ae 
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I. INTRODUCTION 
A. Motivation and background 

Research on ordered condensed matter systems subject 
to bulk random heterogeneities has been an active 
field. Considerable progress has been made, providing 
a better understanding of real materials, where quench 
disorder is always present [H-@- Recently, attention has 
turned to systems where the heterogeneity is confined 
to a surface, e.g. nematic [9l-fll| and smectic liquid 
crystal cells [l2| with dirty substrates. These surface 
disordered systems are of considerable interest and 
exhibit phenomenology qualitatively distinct from their 
bulk disordered counterparts. 

The commonly observed Schlieren texture fl3j is 
a manifestation of such surface pinning in nematic 
cells. Recent studies also include photo-alignment and 
dynamics in self-assembled liquid crystalline monolayers 
[3 d, as well as memory effects and multistability 
in the alignment of nematic cells with a heterogeneous 
random anchoring substrate (l6j . The existence of the 
corresponding phenomena in smectic liquid crystals has 
been recently revealed in ferroelectric smectic-C cells in a 
book-shelf geometry [Tt], EH ■ This latter system exhibits 
long-scale smectic layer distortions, driven by collective 
random surface pinning, and awaits a detailed theoretical 
description. 

A schematic of such a smectic liquid crystal cell is 
illustrated in Fig. [TJ with a "dirty" front substrate 
imposing two types of surface disorders: surface random 
orientational pinning of local nematic order and surface 
random positional pinning of smectic layers. Generically, 
such surface pinning leads to elastic and plastic smec- 
tic disordering, latter characterized by proliferation of 
topological defects, e.g., dislocations. In this paper, we 
focus on the simpler limit of weak disorder, where the 
topological defects are either not present at long scales 
or arc sufficiently dilute, that they can be neglected for 



a range of experimentally relevant length scales. Clearly 
this leaves a rich and challenging regime of strong pinning 
to future studies. 




FIG. 1. A schematic of a half-infinite smectic liquid crystal 
cell in the presence of surface random pinning, with the 
region enclosed by dashed line indicating a domain of area 
£x x £ z (£ z ~ (,x/ a ) within which the perturbative treatment 
of surface disorder is valid. The smectic long range order 
is destroyed by arbitrarily weak surface disorder beyond the 
domain size. 



B. Summary of the results 

We explore surface pinned smectic cell, focussing on 
the bookshelf geometry. As we will show in Sec. IVIIII 
in the homcotropic geometry, the effects of random 
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substrate reduce to that of a random field xy model in 
(d— 1) dimensions, and are significantly less strong than 
in this more common bookshelf geometry. 

Using a harmonic elastic description of smectic liquid 
crystal order with surface quenched disorder, which 
characterizes a thick smectic cell subject to surface 
heterogeneities, we studied the distortion of the smectic 
layers, characterized by a Larkin (Imry-Ma)-like 0, HTJ] 
length scales given by £ x , £ z ~ C^/A that relates the 
distortion along the layer normal (z) and parallel to the 
layer (x) from either random orientational or random 
positional pinning, A = y K/B is the standard smectic 
liquid crystal length, typically comparable to the layer 
thickness a [21| . In the regime of dominant surface 
orientational pinning we find (with Af the associated 
mean-squared pinning strength) 



B 2 X 3 a 2 



(i) 



For the dominant surface positional disorder (of mean- 
squared strength A„) we instead find 



3c- 



1/3 



(2) 



where c = ~ 34.6. In the presence of both types 
of surface disorders, the domain size is determined by 
a combination of these two lengths as discussed later in 
Sec. EVE] 

Adapting the standard Imry-Ma-Larkin analysis to the 
problem at hand, we find that long range smectic order on 
the heterogeneous substrate is unstable against surface 
random orientational pinning in a d dimensional system 
with d < d, , where 



d f — 4 

u lc — 



(3) 



and for the surface random positional pinning in a d 
dimensional system with d < d^ c , where 



dl' c = 6. 



(4) 



Thus, on sufficiently long scales long-range smectic order 
in a three-dimensional cell is unstable to arbitrarily weak 
random pinning of statistically isotropic (nonrubbed) 
substrate. Below these lower critical dimensions ar- 
bitrarily weak pinning destabilizes long range smectic 
order. Within the resulting finite smectic domains 
the correlations of smectic layer distortions grow as a 
power-law of in-plane separation and decay exponentially 
into the bulk of the cell, with the characteristic length at 
depth y > £ x set by y itself. 

We e mploye d the functional renormalization group 
(FRG) |22| - [26l | to assess the effectively strong nonlinear 
pinning physics on scales beyond the smectic domain 
size. The most interesting and potentially experimentally 
relevant prediction of this analysis is a three-dimensional 
(3D) Cardy-Ostlund-like (CO) [27| phase transition at a 



temperature T g , from a weakly disordered smectic for 
T > T g (where on long scales the surface positional 
pinning is averaged away by thermal fluctuations), to 
a low-temperature disorder-dominated smectic-glass for 
T < T g . In the latter phase, the positional disorder 
flows to a CO fixed line at which the correlations of 
layer distortions asymptotically arc characterized by the 
orientational pinning alone, but with an effective strength 
additivcly enhanced with (T g — T) 2 . The high- and 
low-temperature phases are distinguished by (among 
other features) the effective temperature-dependent ori- 
entational pinning strength 



6 f (T) = A f (T)/A f , 



given by 



Sf(T) 




for T>T g 
for T <T n 



(5) 



(6) 



The long-scale smectic layers correlations on the sub- 
strate (at y = 0) are characterized by C(x,z) = 
([uo(x, z) — uq(0, 0)] 2 ), the full behavior of which is given 
in Eq. (|103[) with limits: 



C(x, z) fa a 2 



nx/£, x (T), 



/2^v^76(r), 



for x 3> 
for x <C 



'Xz, 



with the temperature dependent sizes 



(7) 



(8) 
(9) 



These predictions remain valid only on scales shorter 
than the distance between unbound dislocations and the 
scale £f for weak pinning) beyond which the 

nonlinear elasticity may become important. 
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FIG. 2. X-ray scattering peak width (blue curve, as indicated 
by the inset) as a function of temperature. The peak broadens 
below CO transition temperature T g . The phase transition is 
indicated by the flow of random positional disorder strength 
to a nonzero fixed line value A*(T) for T <T g . 
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We also study various experimental features of our pre- 
dictions for the polarized light microscopy. We suggest 
that the highly anisotropic domain (as shown in Fig. Q} 
induced by surface disorder is a plausible explanation for 
the commonly observed stripes in thin smectic cells [l8[ . 
We further argue that the 3D smectic-glass transition 
predicted above (illustrated in Fig. [5]) may have already 
been observed as the aforementioned precipitous X-ray 
peak broadening in cooled smectic liquid crystal cells 
with a random substrate [13, EH > which is expected to 
have a peak width proportional to ^ 1 (T) as shown in 
Fig. [2] Although further systematic detailed studies are 
necessary to test this conjecture, based on the robustness 
of our theoretical prediction, we expect this transition 
to be quite generic in smectic liquid crystal cells with 
heterogeneous substrates. 

The remainder of this paper is organized as fol- 
lows. In Sec. [II] we present a simple model of a 
thick (half-infinite) bookshelf smectic cell with surface 
heterogeneity, and derive the effective surface theory 
by dimensional reduction where bulk modes are inte- 
grated out. in Sec. IIII1 Imry-Ma argument is used 
to estimate the domain size at which the influence of 
random pinnings becomes large. In Sec. IIV1 the surface 
random positional pinning is treated within the Larkin, 
perturbative approximation, that allows us to compute 
exactly the short-scale correlation functions. In Sec. |V] 
we present a functional rcnormalization group analysis 
to treat fully nonlinear positional pinning necessary 
for characterization of long scale correlations. We 
demonstrate that at finite T, this analysis predicts 
a 3D Cardy-Ostlund like phase transition. Possible 
features of observation under polarized microscopy and 
X-ray scattering are analyzed in Sec. [VT] In Sec. I VIII 
the importance of nonlinear elasticity is studied using 
rcnormalization group method to show that it is less 
relevant, thereby justifying the validity of the harmonic 
clastic treatment, and the stability of the oricntational 
order and relevant length scales are also discussed. In 
Sec. IVIIII we derive a theory for a surface disordered 
smectic cell with layers parallel to the substrates and 
show that it mimics a familiar bulk disordered xy model 
in one lower dimension. We conclude in Sec. IIXI with a 
brief summary of the work and future directions. The 
technical details of our calculations are relegated to the 
appendices S [HO El andU 



II. MODEL 

A. Bulk elasticity 

We will focus on a smectic cell in the experimentally 
and theoretically more interesting bookshelf geometry. 
Choosing the coordinate system as illustrated in Fig. [TJ 
so that the smectic layers lie parallel to the (x, z) plane, 
with the average layer normal and director n along the 
z axis. The random substrate is located at y = and 



running perpendicular to the y axis. 

Neglecting elastic nonlinearities [28|, we model a 
half-infinite bookshelf geometry smectic cell with a 
heterogeneous substrate by an energy functional 



H 



bulk 



d A ~ X x 



dy 



f(vi u ) 2 + f(a zU ) 2 



-H„ 



(10) 

in which u(x, y) is the distortion of smectic layer at 
position r = (x, y), and in 3D x = (x,z) spans the 2D 
random substrate, and y is the axis into the bulk of the 
cell. 



B. Surface pinnings 



The surface pinning energy of the substrate is given by 



H„ 



H, 



dp ) 



(ii) 



including uniform pinning of the layers, the coupling of 
the substrate with the directors n and smectic density 
p(x) on the y = substrate, respectively. 

To stabilize the bookshelf cell geometry, a uniform 
component of the pinning is realized on the y = 
substrate via 



(12) 



which dictates the layers to be perpendicular to the 
substrate. In experiment, this is imposed by a uniform 
structureless treatment of the substrate [H, |30| • 

Generically, the interaction between the substrate and 
the liquid crystal also includes the coupling between the 
nematic director h and the local random pinning axis, 
g(x) = (g x , g y , g z ), as well as the coupling between 
the smectic density p(x) and the local random scalar 
potential t/(x). These are determined by the substrate's 
local chemical and physical structure (composition, 
roughness, rubbing, etc...). The nematic director n is 
pinned by the local random pinning axis g(x), with 



Hdn — —Id 



(13) 



We approximated this weak pinning by specializing to 
the smectic state, with the layer normal h taken along 
zVl-Sn 2 +5n and Sn « -V^u @,[I3- With this, H dn 
becomes 



Hdn « J rf d -^| [ 5 2 (x) - 5 2 (x)] (d x u) 2 

+ [g 2 (x) - sj|(x)] (d y u) 2 - 2g x (x)g y (x)(d x u)(d y u) 
+2g x (x.)g z (x)(d x u) + 2g y (-x.)g z (x)(d y u)\. (14) 



For a rubbed substrate used in the experiment [H, [3(| , 
the alignment is highly anisotropic such that \g z \ ^> 
\gx\ ~ \gy\ and pinning of layer normal along z is 
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imposed. Combining with Eq. (|12[) and averaging over 
microscopic scales, we find 



in which the nonlinear variance function R v (uq — u' ) is 
given by 



H y +H dn 
in which 



d d ~h 



(15) 



W = W' + 2[g*(x)-g*(x)] ; 
W Q = 2 F " 



h(x) = -2g x (x)g z (x). 



(16) 



For strong pinning ^-(d y u) 2 , the linear terms in (d y u) 
in Eq. (Tl"4")) can be neglected and d y u taken to be 0. 
Determined by local chemical and physical structure, 
h(x) is expected to be short ranged, with variance 



/ 1 (x)/ l (x') = A / ^- 1 (x-x'), 



(17) 



where 7TT is ensemble average and 5 a d (x) is a (d — 1) 
dimensional short ranged function set by the scale on 
the order of the molecular size a. The precise form of 
6 d T 1 (x) has no qualitative effect on the long scale (longer 
than its range) behavior which is our focus here, and the 
effectiveness of field description of the system already 
averages out details on scales shorter than a. Thus it 
can be replaced by a delta function 6 d ~ 1 (x) without loss 
of generality. 

Similarly, the local random scalar potential U (x) is 
also expected to be short-ranged, and we thereby take it 
to be Gaussian characterized by variance 



I °° p 2 

R v (u -u ) = Au I p 2 a + ^2 y cos { n 9o ["o(x) - uo(x')]} 

\ 71=1 

(22) 

with the first dominant term corresponding to q as 
studied in Sec. |VB] 

Putting these terms together, the surface pinning 
energy is given by 



H 



pin 



x ^f{d x uf - h(x)d x u - V[u (x),x 

(23) 

where we take the random fields h(x) and V[uo(x),x] to 
be characterized by a Gaussian zero-mean distribution 
with variances given in Eqs. (fl7|) and (|2"Tj). and impose 
d y u = on the y = substrate to satisfy the strong 
pinning in Eq. (|15l) . More gcncrically, the random surface 
torque h may also be a function the layer distortion Uq, 
with a periodic variance Aj(m — u' ). However, as we 
will show in Sec. IV Bl such generalization is unnecessary 
at long scales. 



Effective surface model 



U(x)U(X) = Ac/^-^x - x') a A ff ^- l (x - x'), (18) 

with strength Ajj- As studied in the context of bulk 
disorder for smectics in e.g., aerogel @, 0], it couples to 
the smectic density on the y — substrate [with uo(x) = 
m(x, y = 0) the surface layer displacement] 



p(x) = p + J2 Rc{ Pn e m ^ z+U0 ^} (19) 



n=l 



through 



H dp = - j d d - l xp{x)U{x) « - J d d - l xV[u {x),x\. 

(20) 

We also take the random positional pinning potential 
V[uo,x] to be Gaussian with short-ranged variance 



V[ Mo (x),x]^K(x'),x'] 



Pi 



oo 2 

—Re ^e irl! ft)[z-2'+'uo(x)-uo(x')]^ 



xA/^fx-x') 
R v (u - U ' )S d - 1 (x-x'), 



(21) 



Since the surface disorder is confined to the front 
substrate at y = 0, no nonlinearities (within harmonic 
elastic treatment) appear in the bulk (y > 0) of the 
smectic cell. Consequently, it is convenient to exactly 
eliminate the bulk degrees of freedom u(x, y) in favor 
of the layer distortion field on the random substrate, 
uq(x) = u(x.,y = 0). This can be done via a constrained 
path- integral method by integrating out u(x, y) with a 
constraint u(x, y = 0) = uq(x), thereby obtainin g an 
effective (d — l)-dimensional Hamiltonian for uq(x) j3l| . 
Equivalently (for T = properties), we can eliminate 
w(x, y) by solving the Euler-Lagrange equation 



KV 4 Lu-Bd 2 z u = f(x)5(y), 



(24) 



with f(x)6(y) representing the boundary condition that 
imposes uq(x) on the substrate. To this end, we Fourier 
transform u(x,y) over (x, y), obtaining an algebraic 
equation for 



u(q.,q v ) = j d d 1 xdyu(x,y)e~ 



iq-x-iq y y 



(25) 
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whose solution after Fourier transform over q y becomes 
(as shown in Appendix |A"]) 



u(q,y) = u {q x ,q z )e~ 



\/2X 



qi+q 2 z +Xql 



Q 2 



Ml 



V x2( ix + <& - Ml 



2 <li 



^y^^qT^-Xql) 



2A 

uo(qx,qz)</>(<i, y), 



Ml) 



(26a) 
(26b) 



where A = y/K/B is the length that, deep in the 
smectic phase is comparable to the microscopic length 
a (determined by the layer thickness or molecular size) 
and we defined a shape function 0(q, y) that characterizes 
the extension of surface-imposed distortions into the bulk 
(y > 0). In obtaining Eq. ([2"B]) . for strong surface 
pinning (W — > oo) in Eq. (|T5j) . we also imposed a 
constraint d y u = at y = 0. More generally, the relative 
importance of the sine and cosine contributions depends 
on the system and the homogeneous transverse pinning 
strength W. 

After substituting the above solution into (fTCTj) and 
integrating out the y (bulk) degrees of freedom in a 
half infinite system (0 < y < oo), the energy functional 
simplifies to an effective surface energy 



H 



surface 



d d -'q x dg z T q + W Q q x 
(27T) 



d-1 



\u (q x ,q z )\' 



jd-l 



x[/i(x)9 lUo + %,x)], (27) 



confined to the random substrate at y = 0, with 



T q = BV2\y/X>ql + q^ V^q x + q\ + M 2 . (28) 

The resulting effective surface Hamiltonian is convenient 
for studying the surface properties of the system, with the 
bulk properties obtained from Eq. (126f) . Equivalcntly, 
this result can be obtained by solving the Euler-Lagrange 
equation and specifying the boundary conditions, as in 
the nematic liquid crystal cell with surface disorder [ll[ . 

By comparing the induced surface elastic energy T q ~ 
K<lx (iz "C M x ) with surface pinning due to rubbing on 
the substrate, WQq x , we obtain an important length scale 



K 

Wq' 



(29) 



On longer scales, the surface rubbing dominates, and T q 

can be approximated by T q ~ B\/2\q z ■ In this regime 
the effective kernel is given by 



W Q q 2 x 



BV2Xq 



,3/2 



W Q ql, 



with an anisotropic "elasticity" for small q x 
different from the bulk smectic elasticity and T q 
smectic elasticity without a rubbed substrate. 



(30) 

quite 
surface 



However, for a weak pinning Wq this scale can be 
much larger than the sample region and thus Wq can be 
neglected. For the remainder of the paper we will focus 
on this more interesting non-rubbed substrate limit. 



D. Replicated model 

In treating random heterogeneous systems it is often 
convenient to work with an effective translationally- 
invariant field theory. This is possible via the standard 
replica "trick" [Hj], a t the expense of introducing n 
replica fields (with the n — » limit taken at the end 
of the calculation). The disorder-averaged free energy is 



given by F = —Tin Z = —T lim x 



with 



Z n = / [dug 1 ] e 



-Hi 



']/T 



(31) 



where the effective translationally-invariant replicated 
Hamiltonian face [uq] is given by 



H 



(r) 



surface 
n 

E 



dr 2 q x dq z T q 2 



a, p 



— J2 I d d - 2 xdzA f d x \v%(x,z)-u$(x,z) 



a,/3 



2T E / d^xdzRo \v%(x, z) - 4(x, z)\ . (32) 



The advantage of the dimensional reduction in Sec. Ill Cl 
is that formally the problem becomes quite similar 
to the extensively studied bulk random pinning model 
[22II24I [33} in one lower dimension but at the expense of 
a modified long range elasticity encoded in r q , Eq. ([2"5]l. 
As we observe from the coarse-graining procedure, in 
general the random orientational pinning may also be 
field dependent, with A/[uq(x, z) — Uq(x,z)) a generic 
periodic function. We therefore generalized to the form 
in (1321). 



ESTIMATE OF THE FINITE SMECTIC 
DOMAIN SCALES 



III. 



The qualitative features of the response of the smectic 
cell to a heterogeneous substrate can be understood 
through a generalization of the Imry-Ma argument [2(| 
to the surface-pinning problems jol— fill] - In this section, 
we estimate the characteristic length scales beyond which 
pinning energies become significant in distorting smectic 
order. 
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A. Bulk Imry-Ma analysis 



B. Surface Imry-Ma analysis 



For an ordered region of L x x L z on the heterogeneous 
substrate the distortions decay into the bulk within a 
depth of L y ~ L x . Consequently, the elastic energy cost 
of such region scales as 



E e ~ L d - X L Z [K(a/L 2 x ) 2 + B(a/L z ) 2 } 



(33) 



The two elastic contributions balance for L 2 ~ XL Z 
(where A = y K/ ' B) , as an expression of a Virial theorem, 
justified in detail in in Sec. IIVBI With this anisotropic 
scaling the elastic energy reduces to 



E,, 



2 \ 2 rd-3 



Ba X L 



Ba 2 X^/ 2 L{ d - 3 ^ 2 . 



(34) 



In such region, the interaction of the layers with the 
random positional pinning can lower the energy by a 
typical value 

(d-l)/2 



E v ~ V P JN, 



A}/ 2 y/L d x - 2 L e ° 



A l J 2 L d J 2 /V\ 

A l/2 A (d-2)/4 i (d/4) 



(d-l)/2 



(35) 



where V p is the typical random pinning strength with 
zero mean and variance A v w V 2 /^ 1 (£o is the 
pinning correlation length) and N p is the number of 
surface pinning sites. Comparing between E v and E e 
shows that for dimension d < d\ c = 6, arbitrary weak 
random positional pinning always dominates over the 
clastic energy on sufficiently long scales, 



and 



L x >&~ (B 2 /A v y/^- d \ 



L z >6~(5 2 /A,) 2/(6 - d) . 



(36) 



(37) 



Thus, long-range smectic order is always destabilized 
on these long scales by an arbitrarily weak random 
positional pinning. 

Similarly, the random orientational pinning can lower 
the energy through interaction with the layers by a 
typical amount 



E f ~ V h y/Ka/L x 



A l/2 L d/2-l „ A l/2 £ («J-2)/4 



(38) 

Comparing Ef and E e , we find that arbitrarily weak 
random orientational pinning destroys long range smectic 
order for d < df = 4 on sufficiently long scales 



L x >£ x ~ (B 2 /A f ) 



l/(4-d) 



and 



L z > Z, ~ (B 2 /A f ) 



2/(4-d) 



(39) 



(40) 



Reassuringly, as we will see, these length scales, 
identified through the Imry-Ma argument, are consistent 
with the domain sizes given in Eqs. (|6Il I62[) . obtained 
through a detailed field theoretical calculation. 



A complementary (but equivalent to bulk) way of 
estimating the surface-pinning and bulk ordering compe- 
tition is by using the substrate model derived in Sec. Ill Cl 
From the effective surface theory given in Eq. (|2"7| (with 
the pinning of Wq neglected), for an ordered region of 
L x x L z on the heterogeneous surface, the elastic kernel 
r q scales according to 



T q - BX 2 L- 3 - B\ 1/2 L~ 3/2 . 
Therefore the elastic energy cost is given by 



2 2 rd-3 



E e - L a x '^L z B\ z a z L- i - BX z a z L 



(41) 



(42) 



or in terms of L z 

E e ~ L d - 2 L z B\ l ' 2 a 2 L- 3 / 2 ~ BX^ 2 a 2 L^' 2 , (43) 

with the same scaling as Eq. (|3~4"|) . Comparing with the 
scaling of the random pinning energies given in Eqs. (|35[) 
and (|3"8")l . the Imry-Ma argument of effective surface 
Hamiltonian leads to the same critical dimension and 
characteristic lengths as in Sec. IHI Al 

IV. SHORT SCALE "LARKIN" ANALYSIS: 
RANDOM FORCE APPROXIMATION 

A. Random force (linear) approximation 

The importance of surface pinning can be assessed 
by computing smectic layer distortions (uq) (dominated 
by the zero-temperature distortions) within the Larkin 
approximation [l9j, which amounts to a linear random 



force approximation, F(x) = d Uo V[uo(x), x] 



to the 



«o=0 



random potential V A [it (x), x], with inherited Gaussian 
statistics and variance 



F (x)F(x') = A^-^x - x') = -<(0)5 d - x (x - x'). 



In momentum space the pinning energy becomes 
l^q* 



(44) 



pin 



(27T) 



d-l 



iqxh(q x ,qz)—F(q x ,q z ) u {-q x ,-q 



(45) 

Within this approximation, the correlation of the layer 
distortion on the random substrate (y = 0) is 



C 



Larkin 



(q) = 



T A f gl + A v 



r 2 

1 q 



including the thermal contribution 



C 



T.Larki 



T 
1 q 



(46) 



(47) 
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and the contributions from both types of surface disor- 
ders 



r 2 

q 



(48) 



At long scales of interest, the contribution from 
thermal fluctuation is clearly subdominant to the zero 
temperature random pinning driven distortions. We 
therefore focus on this latter disorder dominant contri- 
bution. 



This agrees with our Imry-Ma analysis of the previous 
section and is to be contrasted with the bulk predictions 
of Ref.0, that finds 



(u 2 ) ~ A f 



q 2 d d 1 q x dq z 



(u 2 ) ~ A 



f J i4 + <&) 2 Ji/ Lx 4 

d d ~ 1 q x dq z 



q d x +2 dq x _ L5 _ d 



X i 

(54a) 



(54b) 



1. Power counting 

For various ranges of q x and q z , the elastic kernel T q 
exhibits the following asymptotics 



p f BV2\qz /2 , Ml 
q ~ 1 BV2X 2 q 3 x , Xq 2 x » q z . 



(49) 



with d\^ lk = 5 and 7 for for bulk orientational and 
positional random pinnings, respectively. Thus, similar 
to the surface disordered nematic cell discussed Ref. fiol 
and [n|, nere too the restriction of the pinning to a 
surface reduces its lower critical dimensions down by one. 



B. Domain size 



We note that there is a 2 : 1 ratio in the powers of q x and 
q z , that reflects the anisotropy of the underlying smectic 
state. This anisotropic scaling is crucial in our estimates 
below the critical dimensions for the importance of 
surface pinning. 

The correlation function contribution from random 
orientational pinning is given by 



Larkin 



(BV2\ql /2 + BV2X 2 ql 



(50) 

Thus, the scaling of smectic zero-temperature distortions 
is approximately given by 



f ll/L. 4 
J2 M-2 



q 2 d d 2 q x dq z 4 _ d 



L 



(51a) 



qid a 2 q x dq z 



f hit, <& 



L ( t d)l2 , (51b) 



and leads to the lower critical dimension d{ = 4 for the 
random orientational pinning. 

For the random positional pinning, we have 



A. Larkin 



(?) 



A* _ Ay 

^ f BV2Xqz /2 + BV2\ 2 q% ' 



(52) 



which gives 



( U 3) „ A v I - ^ ~ L 6 ~ d 



(u 2 ) ~ A, 



1/L m Qx 

d d ~ 2 q x dq, 



(53a) 

4 6 - d)/2 , (53b) 



leading to d^ c = 6 as the lower critical dimension for the 
positional pinning. 



The utility of the Larkin approximation is that it 
predicts the range of its own validity, limited to length 
scales where the distortion of u is small. From Eqs. 
(|5"Tj) and ([53")) it is clear that below the respective di c , the 
distortions diverge at long scales, signaling an instability 
of the smectic state and associated breakdown of Larkin 
approximation. 

We identify the substrate extent L x and L z at which 
these smectic distortions grow large to the order of the 
layer spacing a as (the so called) Larkin domain lengths 
£ x and £ z [3, [34|. More specifically, Larkin lengths are 
defined as 



(u%{x,z)) = a 



(55) 



that clearly depends on both surface orientational and 
positional disorder [22T4261 135|. Thus, the values of £ XiX 
depend on the relative strengths of these two types of 
disorder. For L x <C \J XL Z in 3D, a simple analysis (with 
details relegated to Appendix [B| gives 



!(*)} 



A/gg + A. u ^ 1 A„ ^ , 1 A; 



r 2 

q 1 q 



3c B 2 X 3 



B 2 X 3 ' 



A 2 

with c = -2ZL_ « 34.6, leading to the equation for £ x 



(56) 



1 A v 



,3 , 1 A y 



i x = a 2 . 



(57) 



3cB 2 X 3 ^ x cB 2 X 3 
In the opposite limit of L z < L 2 /X we instead get 

7727^Y^ f A/gg + A, _ 1 A, 3/2 1 A f 1/2 
^ W; ~y q T q ~3c5 2 A 3 / 2 2 + cB 2 X^ 2 z ' 

(58) 

that leads to the equation for a £ z 



— A " g 3 / 2 + 1 A/ g 1 / 2 = a 2 (59) 
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As discussed in Appendix [B] above equations are to 
be understood in the scaling sense and not as precise 
quantitative conditions. 

The Larkin lengths £ X)Z in Eqs. ([57)1 and ([5§]l give 
the size of the smectic domains on the random substrate 
within which the elastic energy dominates over the pin- 
ning, which is effectively weak and Larkin approximation 
is valid. 

More generally, in the presence of both types of 
pinning, the smectic domain size is determined by the 
minimum of £| and , calculated in Appendix [B] 



in which 




(60) 



(61) 



and 



S.r 



3c 



B 2 A 3 / 2 a 2 



2/3 



3c^ 



1/3 



(62) 



The highly anisotropic domain constructed by £ x and 
£ z is illustrated in Fig. [1] 




FIG. 3. The correlation functions on the heterogeneous sub- 
strate along and perpendicular to the layers, with dominant 
orientational disorder valid for (x, z) <C On longer scales 
nonlinear nature of pinning needs to be taken into account, 
that we do via RG and matching methods in Sec. FVl 



Surface correlation at short scales 



As shown in Appendix [Cl with Eq. (|46| the C(x, z) 



([uq(x, z) — «o(0, 0)] 2 ) can be easily calculated within the 
Larkin regime (x <C £ x an d z ^ £z)j an d along the 
layer (x) and layer normal (z) direction the orientational 
pinning contribution is given by 



C f (x,0) 
C f (0,z) 



na x 



2ira 



2na 2 y/Xz'/Zl, 



(63) 



where £| z are given in Eq. f|6T|) . At short scales x <C £x 
and z <C £z the behaviors of these correlation functions 
are plotted in Fig. |3l The contribution from random 
positional pinning is 



C v (x,0) 
C v (0,z) 



3a 



87m 2 



(#) ' 



(64) 



in which £j! z are given in Eq. (j6"2")l . At short scales, x <§; 
£ x and 2 -C £ z , these correlation functions are plotted in 
Fig. HI 

Thus, within finite smectic domains (x, z) <C £e, 2 , on 
the random substrate we predict anisotropic power law 
correlations given above. 



D. Estimate of distortions into the bulk 

The nature of the distortions in the bulk of the cell, 
y > 0, can be deduced from the u(q, y) in (|2"6"|) and the 
distortions iio(q) on the heterogeneous surface computed 
above. The main feature of the w(q, y) is that the 
amplitude of distortions at wavevector q = (q x ,<lz) 
exponentially decays into the bulk with length 



£ q = l/£>(q) 



2A 



1/2 



vW+gf+Ag 2 



(65) 



where D(q) = ^ y/^<A+ <ll + Ml 

For surface smectic order limited by domain size 
tenCz) the characteristic substrate surface wavevector 
qL is set by (£~ 1 ,Cj 1 ) an d leads to 

^pecssr 1 )] -1 

1/2 

(66) 



2A 



in which £ x rj \/Af^ is used. Thus we predict that for y > 
£ y (ignoring the more subtle distortions on scales beyond 
the Larkin scale, see Sec. |V| smectic surface distortion 
will anneal away into the bulk. 
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C v (x,0) 

2a 2 - 




FIG. 4. The correlation function on the heterogeneous 
substrate along and perpendicular to smectic layers, with 
dominant positional pinning, valid only inside Larkin do- 
mains, (x,z) <C £e,z. On longer scales an RG analysis is 
necessary. 



This allows us to define a y-dependent smectic domain 
size via Larkin criterion on 



(u 2 (y)) 



(2tt) 2 

xdt 
(2tt)2 



a 2^(y) e -2D[c 1 (y),er 1 (y)]y ; 



(67) 



where we assume dominant orientational pinning (posi- 
tional pinning can be similarly analyzed). 

As usual £ x (y) is set by (u 2 {y)) ~ a 2 , thus we have 

Uy)-^ 2y '^ y \ (68) 

A self-consistent solution for this equation is simply 

Cx(y)~max(^,2y), (69) 

showing near healing of smectic order (but see Sec. [Vj) at 
a distance £ x into the bulk. 



random force (Larkin) approximation is inadequate, and 
the effects of the surface disorders must be treated 
nonperturbatively. As with bulk disorder problems, 
this can be done systematically using an FRG analysis 

MM. 



A. Renormalization group analysis 

We employ the standard momentum-shell RG trans- 
formation [36| by separating the layer distortion field 
into long and short scale contributions according to 
Uq(x,z) = Uq < (x,z) + Uq > (x,z) and perturbatively in 
the surface disorder integrate out the high wavevector 
fields Uq > (x, z) that take support in an infinitesimal shell 
A/b < q x < A = 1/a, with b = e se . We follow this with 
a rescaling of lengths and of the long wavelength part of 
the field in real space: 



x 
z 



bx', 
b w z\ 



u^Jbx' ,b w z') = b^u^{x',z'), 



and in momentum space 



= b~ 
= b~ 



(70) 
(71) 
(72) 



(73) 
(74) 



(b-%,b-"£) = b d -^+t>u«(qi x , q i z ), (75) 



so as to restore the UV cutoff back to A = 1/a. Because 
a smectic liquid crystal is periodic under translations by 
a multiple of the layer spacing, it is convenient to choose 
the arbitrary field dimension = 00 and take w = 2 
so that the associated smectic period and length A = 
\J K/B are not rescaled under the RG transformation. 

For simplicity and clarity of presentation we focus 
on an infinitely thick cell with the Hamiltonian given 
in Eq. ([3"2"f. where both disorder strengths are general 
functions of Uq(x, z), and employ the functional renor- 
malization group method. The above rescaling leads to 
zcroth order RG flows of the effective couplings, that for 
a thick cell are given by 



r' q (6) = & d r q , 

B'{b) = b d - 3 B, 
A' f (u,b) = b d - 2 A f (u), 

R' v (u,b) = b d R v (u), 



(76) 
(77) 
(78) 
(79) 



V. PHYSICS ON SCALES BEYOND SMECTIC 
DOMAINS 

On length scales longer than the crossover scales £ x z , 
the distortions of uq grow into a nonlinear regime, where 



where T is kept fixed, and Ay (it) is a periodic function 
of the field u. 

The Hamiltonian can be separated into three parts 

H^ rface = H a + H Af +H Av , (80) 
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where H Af and -Ha„ are given by 



Ha, = — Y 



a. /J 



id- 2 



tio(x,z) - Tigfoz) 



x <9 S 



(81) 



OL.fi 



(82) 



and H have already been defined above, Eq. ([32]) . 
We limit the FRG analysis to one-loop order, 



J 



performing the momentum shell integration over the 
high-wave-vector components Uq> pcrturbatively in the 
nonlinearities A/[itg (z, z) — Uq(x,z)] and A v [uq(x, z) — 
Uq(x,z)]. The change in the Hamiltonian due to this 



coarse-graining is given by 

SH^ = (H Af ) > + (H Av ) > - 



2T 



(83) 



Relegating the detailed calculations to Appendix [Dl 
this coarse-graining procedure leads to functional RG 
flow equations for A/ (it) and R v (u) 



d t A f {u,l) = (d- 2)A f (u) + V A'}(u) - ^R'l{u)R:\u) - A 4 A f (u)A'}(u) 

~1 



-M [A/WA» - M«)A?(0) - A,(0)A'» 
+A [A»i?» - A'/(0)i?» - A'»<(0)] , 



-A'^A'^-A^A^O) 



flfe-R,, («, t) = dR v (u) - CA / (u) + r)R'l («) + A 



-<(u)<(u)-<( U )<(0) 



+ A 3 



-A/(u)A/(u) 



(84) 



■A/(«)A/(0) - A 5 [A f (u)R"{u) - A/fOK'H - A/(u)<(0)] , 



(85) 



2lT3X> ^ = s^lkpAS and other coefficients 



in which 77 = 

are non- universal constants given in Appendix [Dj 



B. Single harmonic form of positional pinning 

In general, the fully nonlinear functions Af(u) and 
R v (u) in ([32]) need to be treated under the functional 
renormalization group coarse-graining 0, [22T(24| . How- 
ever, at finite temperature and in 3D, it is clear from 
Eqs. ([84]) and ([85]) that at long scales, the pinning is 
dominated by the field independent Aj and the lowest 
harmonic of the positional disorder Q , 



R v {u) = A v cos {q u)/q^ 



where qo = 2ir/a is the smectic wave number. Because 
at finite T, 77 is nonzero, all higher, n > 1 harmonics 
of R v (u) in Eq. ([2"2"]) are less relevant with eigenvalues 
A n = 3 — r/n 2 qQ < Ai and can thus be neglected. 

With these considerable simplifications, the flow of the 
surface disorder strengths in Eqs. ([84")) and ([8"5"]) reduce 
to 



d e A f = A f + ^q 2 Al 

d e A v = (3 - ml) A v - ,4 g 2 A 2 . 



(87) 
(88) 



The dimcnsionless coupling of the random positional 
pinning, A v = Aq$A v = 67r 2 (a/^) 3 which is dimcn- 
sionless measure of the random positional pinning, flows 



according to 



d e A v 



3 - Tg ° ^ A - A 2 



(89) 



Taking an ansatz A v (£) = f(£)e( 3 adSW*, the solution to 
this flow equation can be readily obtained: 



A v (£) 



Til 
2uB\ 



A„(0)e( 3 -^ 



2ttBA 



A„(0) 



e( 3 2irBA )^ — 1 



(90) 



Thus we find that A^ decays to the A* = fixed line for 



(86) T > 6nBX/ql = T g , and flows to fixed line 



A* 



2nB\ J 



(91) 



for T < T g , as illustrated in Fig. [5] 

We therefore predict that a 3D surface disordered 
smectic cell exhibits a Cardy-Ostlund-like [27J phase 
transition at T ff = QirBX/q 2 , between a high temperature 
thermal smectic phase, in which at long scales the 
surface positional pinning is averaged away by thermal 
fluctuations and a low temperature random-pinned 
smectic-glass phase controlled by the nontrivial A*(T) 
fixed line. This surface pinned state (and the associated 
transition) is quite similar to the super-rough pha se of 
a crystal surface grown on a random substrate [371 ] . the 
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N 



V V V 



FIG. 5. RG flow of A„ in d = 3: for T < T a = 6%B\/q%, it 
flows to a fixed line A* = (3 — ^bx ) ant ^ ^ or ^ > T s it flows 
to disorder-free fixed line A^ = 0. This corresponds to a glass 
transition at T g from a high temperature thermal smectic 
phase to a lower temperature randomly-pinned smectic-glass. 



1+1 vortex glass phase of flux- line vortices (confined to a 
plane) in type II superconductors [H|,[38||, and 3D smectic 
liquid crystalpinned by a random porous environment of, 
e.g., aerogel Q- 

The dimcnsionless coupling of the surface random 
orientational pinning, Af = Aq^Af = 87r 4 a/£|, flows 
according to 



diAf = A/ + tt 2 A 2 v . 



(92) 



From this equation, we learn that even if the bare value 
of Af is zero, under coarse-graining, at long scales the 
random field pinning, A,,, generates a nonzero random 
orientational pinning Af(£). Since for T > T g , A v —> 0, 
the flow of Af reduces to d(Af = Af, with a trivial 
solution Af{i) = Afe e . It is clear that the random 
orientational pinning is strongly relevant even in the high 
temperature phase, T > T g , latter distinguished by the 
irrelevance of the positional pinning, A„ — > 0. 

In contrast, in the low temperature phase, T < T g , 
the dimensionless strength of surface positional disorder 
approaches a non-zero value, A„ — > A* > 0. Replacing 
A„ in Eq. (|92[) by its fixed point value we obtain 



A f (£)= A f + ir\A*J 



^ 2 (a;) 



(93) 



Similar to the bulk disorder result Q , as £ — > oo this has 
the same asymptotic behavior as in the high temperature 
phase, Afoe e , but with the enhanced nonuniversal 
amplitude A fg acquiring an additional contribution that 
scales with (T g — T) 2 for T < T g . 



C. Matching analysis 

As discussed earlier, correlation function on scales 
shorter than the Larkin domain size can be computed 



within the random force approximation in Sec. IIVI 
However, at scales longer than the domain size effective 
pinning grows large (compared to the elastic energy) and 
can therefore no longer be treated perturbatively within 
the random force model. Nevertheless, we can utilize the 
RG and matching method to effectively overcome this 
difficulty. It allows us to establish a relation between a 
correlation function at long scales (beyond the domain 
size), which is impossible to calculate directly due to the 
aforementioned infra-red divergences, to this correlation 
function at short scales (below the domain size), which 
can be easily calculated in a controlled perturbation 
theory @, OS H, M M M 53 of Sec. El 

For T < T g , the surface positional disorder is relevant. 
However, near the Cardy-Ostlund transition T g , (3 — 
VQo) ^ 1? an d the fixed point value of dimensionless 
positional pinning is small (A* <C 1). Making use of 
the matching method, we establish a relation between 
replicated correlation functions at short and long scales: 

G a p (^(0)^(0)^(0)] 
Uo(q x ,q z )u^(q' x ,q' z ) 



<5 2 (q + q') 

= e^ l C aP [eWWqW, A } {£), A v (£)] .(94) 

We then choose the value of £ = ln(A/g x ) so that the 
rescaled momentum q x is at the cutoff A and therefore 
the correlation function on the right hand side of the 
above equation can be safely evaluated perturbatively. 
This gives 

C af) [q x ,q z ;T^0),A f {Q),A v {0)] 
= C af} [K,{k/q x y qz ;T*^t),Af{t),A v {t)\ , 

(95) 

and for £ — > oo corresponds to long scales with q x — > 0, 
where £* is big enough so that the disorder strengths on 
the right hand side can be replaced by their fixed values. 

Safely calculating the correlation function at the cutoff 
scale on the right hand side using the Larkin analysis of 
Sec. IIV Al we obtain 



C aP [A, (A/q x yq z ; T* (£*), Af(t), A v (t )] 
T5 aP A f (t)A 2 + A v (t) 



r q o?*) 



Tl(£*) 



(96) 



where the £* dependent parameters the recursion rela- 
tions Eq. (|77H79"|> at £* = \n(A/q x ): 



r q (0 

Af(t) 
A„(0 



A 

q x 

A/(r) _ a 

A{£*)qt q x 
A v {t) = Al_ 
A{£*)ql Aq 2 



^ + A4 {k:f 



(97) 
(98) 
(99) 
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Here A(£) = A is a constant given by Eq. (ID13I) . In 
the last two equations we have taken q x small enough 
and therefore I* large enough so that the A v (£*) has 
flowed to its fixed point value given in Sec. IV Bl and the 
low temperature solution of Af{£) is used, with £* large 
enough that the second term in Eq. (|93|) can be ignored. 
Combining these ingredients we therefore find 



Cap [q x , q 
T5 a/3 



z. r q , Af, a v 



S,T-(I 



a: 



2 



1 A * 
A 3 A^,v 



V 2 
1 q 



-, (100) 



in which T q is the bare long-range surface elasticity kernel 
given in ([25)1 . Af is the bare strength of orientational 
disorder strength, and A* is the fixed value given in 

(p?Tj) . It is easy to see that the %A* term arising from 
the random positional pinning is negligible as q x — > 
compared with the contribution of the surface orienta- 
tional pinning. Thus, the smcctic phonon correlation 
functions in the high and low temperature phases are 
only subtly distinguished by the enhanced amplitude of 
surface orientational pinning from A/ to 



A/(T) = A / 



8n 2 a 



X 



(101) 



A/(T) - A f 




FIG. 6. The exact correction to the long scale dimensionless 
strength of the orientational pinning as a function of 
temperature, obtained by numerically integrating Eqs. (|90l 
I92|l for different "bare" values of the positional disorder: 
A„ = 0.1 (blue solid curve) and A„ = 0.2 (purple solid curve). 
As a comparison, dashed curve is the approximation based on 
the fixed point value of A„ used in Eqs. (T5I [til 153)) . 

The above matching calculation gives the long-scale 
momentum-space surface correlation function 



C(q) 



*fSf(T)ql 
^ q 



(102) 



in which 5 f (T) = A f (T)/A f , as given in Eq. ©, 
distinguishes the high- (T > T g ) and low- (T < T g ) 
temperature phases. 



Moreover, the correction to the orientational disorder 
at long scales, Af(T) — Af, that we approximated in 
Eq. (|92|) by the fixed point value of the positional disorder 
which is zero for T > T g and quadratic in (T g —T) for T < 
T g , is quantitatively inaccurate. The exact solution of the 
flow equation for A / given in Eq. (ftj2")l with A v (£) in ([501) 
leads to a long scale correction of Af in the correlation 

function that depends on the "bare" value of A^ and is 
non-zero for T > T g , as shown in Fig.[B] This difference 
does not change the qualitative behavior of the system, 
e.g. experimental feature of the X-ray scattering peaks 
described in Sec. I VI Bl However, it may be important for 
quantitative comparison with experiments. 



D. Surface correlation function at long scales 




FIG. 7. The value of dimensionless function f(t) with t = 
Xz/x 2 , which has limiting behaviors given in Eq. (|105|) . 



As discussed in Appendix [El the full behavior of the 
real-space long scale correlations can be summarized as 



C(M) = Mi/(t), 



(103) 



4tt 2 £ 2 A 3 

in which t = ^* and f(t) (plotted in Fig. [7]) is given by 

1 

115 1 



/(t) = 7r(7r-2) + 7r 
3tT 



2 + 2e J* spirt — 7rcrfc 



V2 
3Vt 



4' 2' 4' 64i 2 



'l) r (lj lF2 \ 



13 7 1 



(104) 



where iF2(ai; b±, &2; x) is the generalized hypergeometric 
function. The limiting behavior of f(t) is given by 



/(*) 



7r(7r - 2) - iirt 2 , fort<l, 
2MVt, for t > 1. 



(105) 



In the (x, z) plane, the above limits correspond to 
regions near the x and z axes, where simple asymptotic 
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C(x,z) 




C(x, 0) 



FIG. 8. The long-scale surface phonon correlation function 
C(x,z), exhibiting limiting behavior C(x,z) ~ x and *J~z for 
small and large value of t = Xz/x 2 , respectively. 



form of the correlation function C(x,z) is available, as 
shown in Fig. [5J In real space along x, it is given by 



C(x,0) 



i (tt-2)A / (T) 



Att 2 B 2 X 3 2 
= Tra 2 5f(T)x/£l = na 2 x/^(T), (106) 

where ^(T) = £f/5 f {T). Along z we instead find 



C(0,z)^2a 2 5 f (T)^J^ 



(107) 



where £ Z (T) = £[/5 2 (T). Thus at long scales the 
correlation function asymptotically scales identically to 
that in the Lar kin's perturbative regime, (|63[) . but with 
an amplitude enhanced below T g , as illustrated in Fig. [9] 
Within the harmonic elastic theory (i.e., neglecting 
clastic nonlinearities and dislocations) we expect above 
predictions to be asymptotically exact. We explore the 
role of elastic nonlinearities in Sec. rVTI Al where we show 
that they indeed are expected to modify our predictions 
on scales longer than the nonlinear length scale £jvl but 
leave the challenging problem of their study to future 
research. 



VI. EXPERIMENTAL PREDICTIONS 

We next explore the experimental signatures of our 
predictions for the surface pinned smectic cell, focussing 
on the polarized light microscopy and X-ray scattering. 

A. Cross-polarized light microscopy 

In smectic cells with weak random substrate pinning, 
the molecular orientation is expected to vary slowly 



T<r P 




C(0, z) 



T<T, 




FIG. 9. The two-point smectic phonon correlation function 
along x and z axes on the random substrate at y — 0. 
For temperature T > T g the correlation is the same as the 
short scale result, (|63l) . while below T g it is modified by the 
enhanced amplitude Sf(T), ([5l [6j> . 



in the bulk, and under standard conditions we expect 
the Mauguin limit be valid with polarization of light 
following local optic axis. Under such conditions, 
a convenient geometry to probe the substrate-driven 
smectic distortion is that of crossed polarizer-analyzer 
pair, with the polarizer aligned with a well rubbed back 
substrate. In this geometry, the transmission vanishes in 
the absence of smectic distortions. Thus, transmitted 
light directly images the optic axis distortions on the 
front heterogeneous substrate [HI EH- This provides 
a sensitive way to investigate statistics of the surface 
smectic layer distortions d x uo(x, z) in a smectic cell 
weakly pinned by a dirty substrate. 

Here, we only consider smectic order distortions on 
the short scales (within the Larkin regime) of thick 
smectic cells where the effects are most pronounced. 
The subtle long scale asymptotic behavior predicted here 
is unlikely to be accessible in current experiments, as 
it requires currently unavailable detailed quantitative 
analysis. Short-scale treatment below can also be readily 
generalized to long scales. We take the back substrate 
to be rubbed and therefore with the optic (nematic) 
axis well aligned along the z axis, while near the 
front y = substrate the smectic order is subject to 
random oricntational and positional pinning of the dirty 
substrate. 

As alluded to above, in this weakly pinned Mauguin 
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limit optical transmission is determined by the azimuthal 
distortions of the optic axis (layer normal) orientation 
4>o(x, z) = <j>(x, y, z)\ y= o on the random substrate. In this 
geometry with the polarizer at the back substrate along 
the z axis, and analyzer at the front substrate along the 
x axis, the transmitted optical field is given by 



E out (x, z) = E in sm4>o(x, z) 



du Q (x,z) 
'' dx 



This gives polarized light microscopy intensity 
Iout(x,z) sa I in \d x u (x, z)\ 2 . 



(108) 



(109) 



Thus the average transmitted intensity is proportional to 
the variance of the smectic layer normal fluctuations on 
the heterogeneous substrate, 



-y- ss \d x u (x,z)\ 2 . 

J in. 



(110) 



From the local intensity we can also obtain a correlation 
function of the transmitted intensity, related to spatial 
smectic phonon correlations: 

C\x-x! ,z-z') 



lf n \d x u Q (x, z)\ 2 \d x >u (x',z')\ 2 



1 



d x u Q (x, z)d x >u Q (x',z') 
\d x u (x, z)\ 2 



(111) 



where we utilized Gaussian approximation valid at short 
scales. 



1. Surface random orientational pinning 

For a dirty substrate where the orientational pinning 
dominates over positional disorder, the transmitted 
intensity averaged over the Larkin domain can be 
straightforwardly calculated and is given by 



— — = \d x u (x, z)\ 2 

1 in. 



. dq x dq z 4 J_ 
(2tt) 2 9a T2 



(tt - 2) A/ 
4tt 2 B 2 X 3 



l/a 



dq x 



A 



/ 



cB 2 X 3 a 



(112) 



The smectic bulk modulus dependence can in principle 
be probed through its well-understood temperature 
dependence (§1 l4lj. 

In the regime of dominant orientational disorder, 
the transmission intensity-intensity correlation function 



along x axis, C 1 (x, 0) can be straightforwardly computed 
by noticing that 



d x u (x + x',z)d X 'U (x',z) [I 

= a / cos (q x x)dq x 

\d x Uo[x,z)\ z Jo 



and thus we find 



C I (x,0)=I 2 out {l + 2 



sin (x/a) 
x/a 



sin (x/a) 



x/a 



(113) 



(114) 



Similarly, the light intensity-intensity correlation along 
z axis is given by 



C I (0,z) = I 2 out { 1 + 2 



Iz \ \ ira 



(115) 



where Cp(z) = f Q cos (wt 2 /2)dt is the Fresnel integral. 

These correlation functions probe surface orientational 
order, that for weak pinning is long range ordered, with 
structure only on the scale of the layer spacing a. 



2. Surface random positional pinning 

For cells dominated by positional pinning the average 
transmission intensity is given by 



\d x Uo(x,z)\ 2 

(tt - 2)A, f 1/a dq^ 
4ar*BW J 1/(x q 2 



(tt - 2)A„ 
4tt 2 B 2 X 3 



(116) 



where we cutoff the average by the smectic domain size, 
£ X ~Q, given in Eq. (g2J. 

The intensity-intensity correlation function along the 
x axis calculated from (|111|) is given by: 



C I (x,0)*I out {l + 2 



1 + 2 1 



1 f 1 ^ cos (q x x) 

J- / 2 d0 * 

Sa; Jl/£ x 1x 

2 



2& 



(117) 



where we have extended the upper bound of integral to 
infinity 

1 f 1/a cos(q x x) 

— dq x 



1/6* 

+cos( — ) + — Si(— ) 



i 



TTX 

2&' 



(118) 
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valid for < a -C x <C £ x an d where Si(z) = sl " dt 
is the sine integral function. 

Similarly, at constant x along the layer normal the 
intensity correlation function is given by 



From this we obtain: 



S(k)<xJ2pl / dre- %i - ? e tq " z (e l ^ u ^- u ^) (121) 



C T (0,z) 



1 



l/u 



cos (q z z) 



1 + 2 1 



3-4 



2' 



3/2 



7TZ 
i7 



(119) 



valid for < a < z < ( z . The intensity-intensity 
correlations decay to their asymptotic (transmission) 
value on the scale of smcctic domains £ XlZ . 

In fact, we expect the intensity- intensity correlations 
predicted by Eqs. (fl~T4|) , (fTT5]l . (fTT7|) and (jlTO]) to be 
observable in crossed polarized light microscopy on a 
thick smectic cell with a single dirty substrate. 



B. X-ray scattering and characteristic lengths 

Experimental studies of a smectic cell [l7], [HI observed 
that the smectic X-ray scattering peak broadens signifi- 
cantly as temperature lowers towards the A-C transition. 
In this section, we analyze the X-ray scattering utilizing 
the smectic cell model and its analysis from earlier 
sections. 

Following standard treatment of X-ray scattering of 
smectic liquid crystal [4214441 ] . for a wave vector k z 
near k z = ±qo = ±2ir/a, the scattering intensity is 
proportional to the structure function S(k) which is 
the Fourier transform of the density-density correlation 
function. The smectic density p{r) periodic along z can 
be represented by 



where q n = nqy = 2m: /a 



iq„[z+u(r)] 



(120) 



Focussing on the lowest order n = 1 peak at q , and using 
Gaussian approximation for the statistics of u distortions 
(expected to be valid on short scales) we find: 

S{k) oc / dre- i{k *- qo)z - i%± e" 4 (["(^-"(o)] 2 ) . (122) 



For smectic distortions confined to the vicinity of the 
random substrate 1451. the structure function reduces to 



S{k X i kz) 



dxdze- ik * z - ik *- x e—2- G( - x < z) 



(123) 



in which k z — k z — q a . 

A detailed analysis of the peak requires the full form 
of the phonon correlation function C(x, z), which is only 
available at long-scale but is somewhat complicated to 
analyze. However, as discussed by Chen and Toner 
[44j the behavior near and away from the Bragg peak 
is respectively controlled by long and short length scale 
correlations. 

For k x < V / ^ 1 k z , the main contribution to the 
structure function comes from x > s/Xz > £ x , thus 
approximation C(x,z) rj 7ra 2 x / '£ X (T) can be used. At 
finite k z we thereby have 



S(k x , k z ) ~ S x (k x ) 



dxe-^e- 2 " 3 */^ 

4^/6 (T) 
kl+A^[UT)}- 2 ' 



(124) 



where the peak width is given by Ak x = 2tt 3 [£ x (T)] 1 oc 

For a/a- 1 k z < k x , similarly, long-scale phonon correla- 
tion can be approximated as C(x, z) w \f2~Tta 2 U z /£ z (T) , 
and at finite k x the structure function 



S(k x , k z ) ~ S z {k z ) 




1-2F C 



2n 2 



} , for KUT) « 1, 



r 



Sill 



2tt 5 



kMT) 



1-2F S 



2n 2 



(125) 



in which Fc(x) and Fg(x) are Fresnel integrals. Thus we find that S z (k z ) peak width is given by Ak z oc 
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S x (k x ) 




FIG. 10. Structure function from Eqs. ([T24]) . (fT25]) and ([127)1 . 

showing X-ray scattering peak profile at finite k x and k z for 
T > T g and T < T 9 , in which the position of the plotted 
profile function is shown on the Bragg peak by red dashed line 
in the inset of each figure. The peaks broaden at temperature 
below T g . 



MTT'cxSjiT). 

Above details of the peak profile are summarized 
in Fig. [10] We thus predict that the X-ray peak 
widths along k x and k z will be broaden with reduced 
temperature below the CO phase transition temperature 



T g . This prediction for the finite smcctic order charac- 
terized by £, x ,z{T), may have already been observed as 
the aforementioned precipitous X-ray peak broadening 
in cooled smectic liquid crystal cells with a random 
substrate [13, EH- More detailed experimental studies 
are necessary to verify this conjecture. 

Away from the Bragg peak, the structure function is 
dominated by short-scale correlation, well approximated 

q 2 2 

by e ~^ ctyX ' zS> « 1 — ^C(x, z). At large k x and k z we 
thus obtain 

S{k x ,k z )™\c A {k x X), (126) 

with limits 

S(k x , 

^ g 2 A f kl + A v \/^ki + ~k 2 z ~Xkl 
~ 2 2B 2 X fcf (A 2 fc* + ifef ) 

i£a£ If for fc x « VA- 1 /^-, 
as illustrated in Fig. [TDJ 

VII. LIMITS OF VALIDITY 
A. Analysis of elastic nonlinearity 

So far we have limited our analysis (see Sec. [TTJ) of 
a smectic cell with a random substrate to a harmonic 
clastic description neglecting the nonlinear elasticity. 
However, these are known to be imp ortant in smectic 
systems in the presence of thermal [28l |46| and quenched 
random @,S3] fluctuations. In this section, we formulate 
the treatment of these nonlinearities and estimate their 
effects in surface disordered smectic cells. 

A complete description that includes elastic nonlinear- 
ities and random surface pinning (latter limited to Larkin 
harmonic description) is captured by the replicated 
energy functional 
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d z u a - i (V ±U Q ) 2 



i n f r 

"j^fE d d ~ 2 xdzdyd(y) A f d x u%(x,z)d x u%(x,z) + A v u%(x,z)u%(x,z) 

a. 13 

= Hj d d - 2 xdz J™ dy f 



ct,/3 



9f^2 I dd 2xdz dyS(y) Afd x ug(x,z)d x u^(x,z) + A v Uq(x,z)u^(x,z) 



-E/ d d - 2 xdzdy^ {d z u a ){V^u a f + Y, J d d - 2 xdzdyj(V±u a ) 4 



Ho + Ht. 



(128) 
(129) 



I 

To assess the influence of -Hi nonlincarities we perform correlation function derived in Sec. IIV 
(d — 1) dimensional coarse-graining at fixed y, pertur- 
batively in Hi, which allows us to utilize the harmonic 



C{q x ,q z ,y) 



(\u(q x ,q z ,y)\ 2 ) 

T A f ql + A v \^ 

7^o a ,p H — 



)tf(q,y), (130) 



where </>(q, y) is given in Eq. (|26|) . The lowest nontrivial 
correction is given by 



SH = -^(HDl = - / d d - 2 xdydzd d - 2 x'dy'dz'[d z u a (x,y,z)} [d z u^x',y',z')] 

x (((Viu a (i, y, z)] [V' ± uP(x', V, ^')]>>) 2 , (131) 



where the correlator of short scale modes is given by [sec (11301) ] 



([d x u"(x,y,z)\ [d x ,uP(x',y',z')}) > 

d d - 2 k x dk z (T5 af} A f kl + A v \ l2 



dk z ( . 



(2tt) 

x e ^( w )^(k,y)# iy '). 



(132) 



Redefining variables X = (x + x')/2, Sx = x — x', Z = 
(z + z')/2 and 5z = z — z', the dominant contribution is 
given by 



SH =\zZ J d d - 2 XdZ J dydy'8B{y,y') [d z u a (X , y , Z)} [d z u a (X, y' , Z)\ , 



(133) 
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where 



SB(y,y') 



B 2 
2T . 

B 2 f d d - 2 k x dk z ( T A f k 2 x + A„ ^ 



d d " 2 SxdSz (([V_LU Q (a;, y, z)] [V' ± u^(x', y>, z>)])> 

fc^ 2 (k,y)0 2 (k,y') 



2T J (2n) d - 2 2tt Vr k 
-B 2 [ £^h$l A f k l+ Av ki<f> 2 (k, y)0 2 (k, y>). 



(27T) 



6,-1 



2tt 



(134) 
(135) 



Fourier transform with respect to X and Z , ()133|) gives 



5H = lJ2 



d d 2 q x dq z 
d-i 



(27T) 

</>(q, y)0(q, y') 



<W<5B(y,y')CK(q)l 



(2 7 r) d - 1 ^vq)cl"o(q)l > 



(136) 



in which the momentum dependent correction to the 
smectic compressional modulus in the effective surface 
theory, SB(q), is given by, 



6B(q) 



dydy'5B(y, y')<f>(q, y)(f>(q., y) 



= -B u 



d d ~ 2 k x dk z Afk% + A v 
(27r) d - 2 "27 



fc^ 2 (q,k), 



with the function 6(q, k) given by 



6(q,k) 



dy(j) 2 (k, y)0(q, y). 



(137) 



(138) 



We thus obtain the effective randomly surface pinned 
smectic model that includes the effect of smectic non- 
linear elasticity. Although formally well defined, so far 
we have not succeeded in systematically treating this 
correction analytically. We leave such analysis to future 
studies. 

Focusing on the leading exponential contribution 
allows us to carry out y and y' integrals, giving 



6B(q) 



-B 



d k^dk 2, 



A f kl + A v fc4 



(27r) d - 1 Tl [2D(k) + £>(q) 



2 n x> 



(139) 

where D(q) = -^y^/A 2 ^ + q 2 + Xq 2 , and thus in the 
q — >• limit it gives 

x r B 2 f d d - 2 k x dk z A f k 2 + A„ , 

SB {q ->0)« SBq ~ -— J I,' ^ ^ 



(140) 



„ A l L 5-d + A JL L 7-d 

B ^ B x 

Equating SB q 2 to the harmonic kernel in Eq. 
allows us to define nonlinear crossover length ^nl,x 
beyond which nonlinear elastic effects must be taken 



into account. Focusing on orientational pinning (that 
we showed is dominant at long scales, see Sec. [Vj we find 



SB q 2 z \ qz= ^ NL z 
A/ 



B\ 2 



5-d-4 



which gives 



NL,x 



(X 2 /A/) 1/(4 " d) , 



(141) 



(142) 



as the length above which the nonlinear elasticity 
dominates |48j . 

Based on this we expect that although nonlinear 
elasticity is relevant below d^ L = 4 for the orientational 
pinning, it is only relevant at length scales much larger 
than the (Larkin) domain size £, x . Thus, to study system 
properties at length scales smaller than the nonlinear 
length £nl,x, it is safe to ignore the nonlinear elasticity. 
Here we employed this approximation and left the full 
analysis of nonlinear elasticity for future study (49j . 



B. Stability of orientational order 

All of the above analysis of smectic (positional) order 
was predicated on the assumption that the long-range 
ncmatic (orientational) order is stable. Below we 
examine this assumption. 

In the limit of dominant orientational pinning, the 
mean-squared distortion of orientational order on the 
heterogeneous substrate at y = is given by 



(|<Sn(x)| 2 ) « (|a^(x)P 



dq x dq. 



z „4 



(^-2)A y 
47r 2 B 2 A 3 



l/a 



dq x 



l/L 



A 



cB 2 X 3 a 
ali{ « 1- 



(143) 



Evidently, orientational order is indeed long-ranged for 
weak orientational pinning. 

In contrast we find that for dominant positional 
pinning, the distortion of nematic director diverges at 
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long scales according to 



(\6n(x)\2}n(\d x u (x,zW) 

_ (ti--2)A„ f 1/a dq x 



4tt 2 5 2 A3 Jl/L q < 
{n-2)A v 



4tt 2 5 2 A 3 



This implies that the orientational order is destroyed at 
scales above 

fo,» = ^- = i(S) 3 /a a »C (145) 

set by the condition on L such that 5n ~ 1. Similarly 
along z, the length scale beyond which the mean squared 
orientational distortion (|<5n(x)| 2 ) grows to order 1 is 



(146) 



Naively this purely harmonic calculation implies that 
the analysis and prediction of previous sections are only 
valid up to this finite nematic domain scale, £o.x an d 
£o,z- However, based on Eq. (|100l) and the supporting 
RG and matching analysis, we note that at long scales 
the positional pinning becomes subdominant to purely 
orientational pinning, whether it is relevant for T < T g 
or irrelevant for T > T g (in the former case just modifying 
the strength of the effective orientational pinning) . Thus, 
we conclude based on this and Eq. (|143|) that the 
orientational order is indeed stable for weak smectic 
surface pinning. 



VIII. SMECTIC CELL WITH HOMEOTROPIC 
ALIGNMENT 

So far our focus has been on a smectic cell in 
the bookshelf geometry, corresponding to molecular 
alignment parallel to the substrate, with the smectic 
layer normal lying in the plane of the substrate. As 
we demonstrated in previous sections, random surface 
pinning in this geometry has significant qualitative effects 
on the smectic order even in the bulk of the cell, 
destabilizing it on random substrate. 

In this section, for completeness, we also analyze 
a cell with a homcotropic alignment, with smectic 
layers running parallel to the substrate. As we show 
below, in this geometry the problem reduces to that 
of surface-pinned xy model. It thereby maps onto the 
previously studied plan ar nematic order with surface 
random pinning fiol. 

In this homcotropic geometry the cell is described by 
the Hamiltonian 



H sm —Id 



d-l. 



dz 



f(VitO a + f(fl,t*) a 



H. 



pin ; 

(147) 



with the layer normal z perpendicular to the substrate. 
Substrate randomness introduces distortions in the pre- 
ferred layer normal 



(Sn x , Sn y , z) I yjl + 8n% + <5n. 2 , (148) 



(144) leading to pinning energy 



H,, 



= / d d - 1 xdz5{z){-W{h- zf 



g(r)] 2 -y(u,r)} 



d d - l xdz5{z) 
V(u,x 
d d ~ l xdz8{z) 



W + —{8nl + 5n 2 y ) - h(x) ■ 5n 



W 



(V_lw) 2 - h(x) • (V ± u) - V(u,x) 

(149) 



Above we made use of the approximate relation between 
gradient of the smectic layers and the nematic director 
deviation <5n = (Sn x , 8n yi 0) ~ V_lu, and in the 
last equation ignored irrelevant constants. The first 
pinning term encodes the average substrate homcotropic 
alignment, [n ■ g(r)] 2 ~ h(x)-(Vj_ii) captures the random 
orientational pinning on the substrate, and the last term 
is the random positional pinning at z = 0. 

As for the bookshelf geometry, it is convenient to 
carry out the dimensional reduction to the surface field 
u(x,y, z)\ z= q = uo(x,y) at z = in a half-infinite cell. 
To this end we solve the corresponding Euler-Lagrangc 
equation subject to the substrate boundary condition, 

KV\_u{x, y, z) - Bd^u(x, y, z) = 0, 

whose solution for a thick cell is given by 

u(q,z) = u (q)e- x « 2z . (150) 

Plugging this back into the bulk Hamiltonian and 
integrating over the z degrees of freedom for z > 0, we 
obtain the effective surface Hamiltonian 



TJ 



surf 



d 2( l To 2| / \|2 

-q K(q)| 



(2tt) 2 2 

J d^xdzSiz) [h(x) • (V ± u) + V(u,x)] 



(151) 



where 



T = VKB + W. (152) 
In the momentum space the random pinning is given by 
d 2 q 



H 



pin 



[iq ■ h(-q)u (q) + /(-q)"o(q)] 



(153) 
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Thus as anticipated in the Introduction in this 
homeotropic geometry the effects of the random sub- 
strate in a smectic cell reduce to that of a well-studied 
random-field 2D xy-model H S3 HI M, HI for the 
displacement of the smectic layer in contact with the 
random substrate. In a weak pinning limit (where effects 
of dislocations and of nonlinear elasticity can be safely 
neglected) it is predicted to lead to weak breakdown 
of long-ranged smectic order with log-squared phonon 
correlations [H, [37], . The random substrate pinning 
is significantly less disruptive in this geometry and can 
be used to test numerous theoretical predictions for the 
2D random field xy-model QUEH- 



IX. CONCLUSIONS 

A. Summary 

Motivated by a number of liquid crystal experiments 
on smectic cells [13, EH], m this paper we studied 
the smectic liquid crystal order perturbed by a ran- 
domly heterogeneous substrate. Demonstrating that 
the resulting random surface pinning comes in two, 
orientational and positional forms, we introduced a 
general model for a heterogeneous smectic cell and 
analyzed the response and stability of its smectic order 
to a random substrate. We thereby demonstrated 
rigorously that smectic long-range order in a thick 
cell (without a rubbed substrate) at long scales is 
unstable to arbitrarily weak surface disorder, replaced 
by random smectic domains whose size and correlations 
on the heterogeneous substrate and in the bulk we have 
computed. We showed that for weak random pinning the 
smectic order is thereby replaced by a three-dimensional 
smectic glass-like state, that undergoes a temperature 
and pinning-driven Cardy-Ostlund-like phase transition 
between a high-temperature weakly-pinned smectic and 
a low-temperature smectic glass, where smectic layers 
in bookshelf geometry are strongly pinned by the dirty 
substrate. 

We computed the statistics of the induced positional 
and orientational distortions within and on scales beyond 
the smectic domains. We also analyzed the transmission 
signal in the polarized light microscopy and the X-ray 
scattering structure function in such a smectic cell. 
We showed that these can be used to experimentally 
probe our predictions of the finite smectic correlation 
length and the statistics of the substrate induced 
smectic textures. Based on this, we suggested that the 
precipitous broadening of the X-ray peak with reduced 
temperature observed by Jones and Clark [13, EH may 
be associated with the predicted smectic glass phase 
transition. More systematic experimental and theoretical 
studies are necessary to explore this further. Finally, we 
also demonstrated that in the geometry where layers run 
parallel to the random substrate, the distortions on the 
substrate are governed by the well-studied 2D random 



field xy-model. 



B. Future directions 

While our current study makes significant progress 
toward understanding of randomly surface-pinned smec- 
tics, it also stimulates a number of interesting open 
questions, that we leave to future research. Probably 
the most challenging of these is the nature of the locally 
smectic state in the presence of strong random surface 
pinning. This regime will exhibit distortions beyond 
the elastic approximation, resulting in the proliferation 
of dislocations that are notoriously difficult to treat 
analytically, requiring sophisticated numerical analysis. 

A more systematic (e.g., RG) treatment of nonlin- 
ear elasticity, that we only briefly touched on via a 
self-consistent scaling analysis, is another challenging 
problem requiring a simultaneous treatment of pinning 
and elastic anharmonicities. 

Throughout this manuscript we have also assumed 
that the average strain of the pinned smectic state 
vanishes. While this is true in equilibrium, experimental 
observations fl7l . E§1 . [52j suggest that the system may 
have difficult time equilibrating as the bulk layer spacing 
changes significantly with temperature, while the layers 
remain pinned at the random substrate. Thus, as the 
system is cooled into the smectic state, we expect a 
significant build up of strain near the random substrate, 
with relaxation rate limited by nucleation and motion 
of edge dislocations. It is clear that accounting for this 
background strain is crucial for a detailed understanding 
of current experiments (l7l [Tp. We are currently 
exploring this interesting effect [49f . 

Finally, our earlier work on nematic cells [l(| [H| an d 
current study of smectic cells [l2j lend themselves to 
generalizations to a broad array of other liquid crystal 
phases, such as for example the cholcsteric and smectic-C 
states. The effect of the random substrate depends 
qualitatively on the nature of the phases and we leave 
their analysis to future studies. 

Wc hope that the present study will stimulate further 
theoretical and experimental work on these and many 
other open questions. 
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Appendix A: Solution of the Euler-Lagrange 
equation 

In this appendix, we study the solution to the 
Euler-Lagrange equation of the smectic liquid crystal 
cells subject to a boundary condition of a random surface 
texture u(x,y = 0,z) = Uo(x, z) or equivalently random 
force f(x, z) as in (f2"4"|). 

KV\u - Bd 2 z u = f(x, z)6(y), 

where we used a surface pinning force 

Hpi n = - J dxdydz6(y)f(x, z)u(x, y, z), (Al) 

to encode the boundary condition u(x,y = 0, z) = 

U (X, Z). 

In Fourier space Eq. (|24|) reduces to 



[K(ql + ql) 2 + Bql] u(q x ,q z ,q y ) = f(q x ,q z ), (A2) 



giving 



u(qx,q z ,q y ) = 

Then, for y > we have 
dq y 



f(qx,q z ) 



K{ql 



Bql 



(A3) 



u{q x ,q z ,y) 



, 2tt K{ql 

f(qx,q z ) 

2ttB 
f{qx,q z ) 



dq y 



Bql 

e iQyV 



2irB\/\ 



dq y 



V(q 2 + q 2 y ) 2 + q 2 

e tq y V 

{\ql + ql) 2 + qV 



(A4) 



with the associated length A = y/ K/B, q y = \f\q y 

and y = y/y/X. As shown in Fig. [TTJ (|A4j) can readily 
evaluated by a contour integration, closing the contour 
in the upper half of the complex plane. Defining q p = 

\J—\q 2 + iq z , the four poles are given by ±q p and ±q*. 
Closing the contour in the upper half plane we find 



i>{qx,qz,y) = 



f{qx,q z ) 

2ttBV\ 

f(qx,q z ) 

2nB\/\ 

f{qx,q z ) r 

2B\/\ ' 

f{qx,q z ) 



dq v 7= 



2m 



(q y - q P )(q y - q* P )(q y + q P )(q y + q P ) 



(qp - q P )(q P + q P )(q P + q* P ) (-<?; - q P )(-q* P - q* v ){-q* P + q P ) 



-J'-lvU 



4iq p Rc(<7 p )Im((jp) Aiq*Re(q p )lm(q p ) 



2By/X 2Re{q p )lm(q p ) 



Re 



=*<?p2/ 



>ll> 



(A5) 



and 



2 \ 2 

lp = ~Mx 



„ i \tt— arctan ( 

= \ql\e L V^ 3 



(A6) 



with \q\\ = \J\ 2 qi + q 2 - Noting that 

I \2 4 i 2N 1 / 4 Mf _ 5 £ 

q p = (A q x + q z ) e L 
Re(g p ) = (A 2 g4 + ^ 2 ) 1/4 s in 



We find 



1 I qz 
— arctan — ^ 

2 KM 2 



Im(q p )=(\ 2 qt + qt) 1/4 > 



Ml, 



1 ( qz 

— arctan — ~ 

2 \Mi 



l j=\] y/X 2 qi+ql + \q 2 x 



V2 



(A7) 



(A8) 



(A9) 
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u(g x ,qz,y) 



f{qx,q z ) 



2BV2X q z ^\ 2 qi + ql 



y 



\ql cos (-l/v/^-^) 



+V VAM + 9f + Agg sin 



A<?< 



(A10) 




+ * <7z 



V 1 ^ 



«<7z 



FIG. 11. The integral (|A4|) represented by contour integral 
closed in the upper-half of the complex plane with poles 
shown. 

On the y = substrate, 

Uo(q x ,q z ) = u(q x ,q z ,y = 0) 



(All) 

Combining this with Eq. (|A10[) allows us to express 
u(q x ,q v ,y) in terms of this substrate deformation 



u(q x , q z , V) = u (q x ,q z )e-^^ X2 ^ 



J ^\ 2 q x + q 2 + \ q 2 

— S1I 

V x2( ii + iz- Ml 



Ml 



- cos 



<2\ 



2 lt + 9z _ Ag 2 



given as Eq. (|26[) in the main context. 



Appendix B: Domain size in the surface disordered 
smectic cell in bookshelf geometry 

In this appendix we present the details of the analysis 
of the Larkin lengths setting the domain size and shape 
within which smectic order is stable and the weak random 
pinning is treated perturbatively. To make analytical 
progress we focus on extreme limits of dominant posi- 
tional pinning A„ or dominant oricntational pinning A f , 
but not both. The calculation here is only precise in the 
scaling sense and is not quantitatively accurate. 



1. Larkin lengths for dominant surface 
orientational pinning 

Focussing on the case of dominant random orien- 
tational disorder and generalizing the analysis to d 
dimensions with (d — 2) dimensional x, we find that the 
smectic phonon variance is given by: 



(ul(x,z)) = 



A 



2B 2 A 



d 



d-2 



qxdq z 2 V X2 qi +qj- Ml 
(27T)"- 1 Qx q 2 (\ 2 q$ + q 2 z ) 



(Bl) 

This integral is divergent at small q and needs cutoffs 
set by the domain size L x and L z . For L x <C \/\L z , 
one can integrate q z out, cut off the IR integration by 
q x ~ I/Ac j an d find 



A 



/ 



jd-2 



q x it 



4ttB 2 \ J (2w) d - 2 X 2 ql 

( 7 r-2)A / s d - 2 _ r 1 ^ q 

7 n/L 



4irB 2 \ 3 {2n) d - 2 



d-3 



dq x 



47rB 2 A 3 



Cd- r < 



1 



-L 



4-d 



4-d 



(B2) 



diverging for d < d\ c = 4, where Sd is the surface area of a 
rf-dimcnsional unit sphere, and Cd = Sd/(2ir) d . Defining 
the Larkin domain size £f as the length scale over which 
u rms is on the order of smectic layer spacing, we find 



if 



4(4-d)7r J5 2 A 3 a 2 



2)C, 



d-2 



l/(4-d) 



(B3) 
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that for the physically relevant case of d=3 is given by 

B 2 \ 3 a 2 



S.J; 



(B4) 



in which c = ^ « 34.6. 



At the lower critical dimension of d = 4, it rms diverges 
logarithmically leading to 



^ = aexp 



4^ (2^) 2 B 2 A 3 a 2 l 



2 5 2 



(B5) 



Similarly for the Larkin length along z, integrating over 
q x first, we find 



{ul(x,z)) 



A f fdq z f d d ~ 2 q x 2 ^X 2 q x + g 2 - Ml 



2B 2 \ J 2tt J (2?r) 



- 2<7 * q 2 {X 2 q% + q 2 ) 



A/ S d _ 2 /"dfe [°° ^JX 2 q x + q 2 z -Xq 2 x d _ x 



2B 2 X(2n) d - 2 J 2ir J q 2 z (X 2 q x + q 2 z ) 

A f S d - 2 r(i-f)r(|)-^ 3 / 2 sec(f ) 



2iTB 2 X d / 2 + 1 (2ir) d - 2 

f 1/a is 
x / dq z qi 



40F 



1 



A, L 2 - d/2 



2B(d) B 2 X d / 2 + 1 2-d/2' 



(B6) 



in which 



i _ i s d - 2 r(i-f)r(f)-7r 3 / 2 sec(f ) 



B(d) 7T (27r) d - 2 
This gives 



40F 



(B7) 



£ 



/ _ 



(4 - d)B(d) 



B 2 \ d ' 2+1 a 2 



2/(4-d) 



(B8) 



which in the physically relevant three dimensions reduces 
to 



B(3) 



B 2 X 5 / 2 a 2 



A 



f 



(B9) 



where B(3) w 37.5 is roughly the same as c. At the lower 
critical dimension we instead find 



£{ = aexp 



2 B 2 X 3 a 2 



[B(A) A f 



(BIO) 



2. Larkin lengths for dominant surface positional 
pinning 

For a smectic cell where positional disorder dominates, 
the variance of smectic distortions is given by 



if 



3 dq z A v ^X 2 q x +q 2 -Xq x 
c 2n2B 2 X q 2 (X 2 q x + q 2 ) 
dq x it - 2 



AitB 2 X J 2ir X 2 q% 
A v 2(tt-2) f 1/a dq^ 

Jl/Lx Qx 



4ttB 2 X 2irX 2 



(tt-2)A„ 3 



(Bll) 



12tt 2 B 2 A 3 

which diverges as L 3 in 3D. 

Consequently, the Larkin length along x is given by 

S 2 A 3 a 2 iV3 



£ r 

Sx 



3c- 



(B12) 



Similarly, the Larkin length along z is instead given by 



A(3) 



ff 2 A 3 / 2 a 2 



2/:! 



(B13) 



where ^4(3) 
as c. 



39.83, roughly the same 



by 



i6r 2 (f)-V2^3 

More generally in d dimensions the variance is given 
A v f d d - 2 q x 7T - 2 



{u 2 {x,z)) 



AttB 2 X J (2ir) d ~ 2 X 2 q% 

(tt - 2)A„ S d ^ 2 fV a qt 3 dq : 



4nB 2 X 3 (2n) d - 2 J 1/Lx q x 

- 47 r£ 2 A 3 ' (BW) 

which is diverging for d < dj = 6, leading to the Larkin 
length along x given by: 



Sa; 



4(6 -o> B 2 X 3 t 



l/(6-d) 



.(7T-2)C d _ 2 A„ 

And similarly along z 

, v _ \{Q-d)A{d) B 2 X d ' 2 a 2 ^/ 



(B15) 



3 A„ 
in which 

i _ i 5 d _ 2 r(i-f)r(^)-7r 3 / 2 csc(f ) 



(B16) 



A{d) tt (2tt) 



d-2 



12V^ 



. (B17) 

At the lower critical dimension d = 6, u rms diverges 
logarithmically leading to 

4tt (2 7 r) 4 J B 2 A 3 a 2 " 



££ = aexp 
£~ = a exp 



TT — 2 S4 At, 

2A(6) B 2 X 3 a 2 ' 
~ AT" 



(B18) 
(B19) 
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Appendix C: Calculation of correlation functions at 
short scales (Larkin regime) 

Here we present the details of the correlation function 
of smectic phonons Uq(x, z) on the dirty substrate, y = 0, 
focussing on the short-scale Larkin regime. In 3D it is 
given by a Fourier transform of Eq. 

C(x,z) 



(Mx,z)- Uo (0,0)] 2 ) 

2 r ^^ A /g g + A _ eiqxX+i 

J 2tt 2tt (r„)- 



where within the Larkin regime the integrations are 
bounded by the intermediate length scales £ XtZ given in 
Sec. IIVBI We treat the two forms of orientational and 
positional pinnings individually, for simplicity assuming 
one or the other (but not both) dominates. 



1. Dominant surface orientational disorder 

With dominant surface orientational disorder the 
correlation function is given by 

C f (x,z) 



(Mx,z)-u o (0,0)] 2 



f dq z dq x A f ql y/ \ 2 q* + q 2 z - Xq 2 x _ iqxX+iq ^ 
J 2n 2ir 2B 2 \ q 2 (X 2 q x +q 2 ) [ >' 

(C2) 

For the dependence along x, we first integrate out q z 
obtaining (for x <C £ x ) 



C/(x,0)«2 
w 2 



(n-2)A f . 

q x [1 - cos {q x x)\dq x 



4n 2 B 2 X* J ml 
(tt — 2) A/ ttx 



Att 2 B 2 X 3 2 
~ ^a 2 x/^ x . 

For the dependence along z, we instead have 

1- cos fez) 



(C3) 



C/(0,z)«2- 



2 J mi 



dq z 



3/2 



(C4) 



where the definition of z is given in Sec. IIVBI 
These correlation functions with dominant orientational 
disorder are plotted in Fig. [31 



2. Dominant surface positional disorder 

With dominant random positional pinning the corre- 
lation function in 3D is given by 

C v (a;, z) 



[u (x,z) - uo(0,0)] 



= f dg z dq x A v \/X 2 q* + q 2 z -Xq 2 _ iq^x+iq^ 
J 2tt 2n 2B 2 \ q 2 z (X 2 q x + q 2 ) [ ' h 

(C5) 

Focussing first on the x dependence at z = 0, we integrate 
out q z and find 



C v (x,0) 

(tT-2)A. u 

~ 4ir 2 B 2 \ 3 

(tt-2)A„ 
4tt 2 £ 2 A 3 



[1 - cos (q x x)} dq x 



l/x 



-4 (ixX) 2 

Qx ~^—dq x 



<l x dq x 



l/x 



(tt-2)A 2 
4tt 2 B 2 A 3 fe ' 



(C6) 



where we simplified the integral by separating it into 
q x < l/x and q x > l/x parts with corresponding ap- 
proximations. Utilizing the condition from Sec. IIV Bl and 
Appendix |B] that ^Tg^a = 3a 2 /(££) 3 , the correlation 
function reduces to a simple form 



C v (x,0) « 3a s 



(C7) 



Equivalently, this integral can be evaluated exactly in 
terms of special functions and then small x <C tj, x limit 
taken. Consistent with the above estimate, we find 



C v (x,0) 

(7r-2)A„ 
~ 47r 2 B 2 A 3 



1 - cos {q x x)]dq x 



+ \q [x 2 2(0 2 ] cos(|) 



4tt 2 B 2 A 3 I 3 



+ -x(C T ) 2 sm( — )- — 
6 12 



3a z 



3 r 



7T-2 Si( — ) 



(C8) 



where Si(z) = J Q 2 sin (i) ^ is the sine integral function. 

To obtain the short-scale z dependence of the correla- 
tion function, we instead first integrate over q x , obtaining 



C v (0,2) 



a„ ier 2 ( 



2/5 ' 



/27T 3 



4tt 2 B 2 A 



2V^A 



5/2 



2 A„ 16r 2 (|)-V2^2^ _3 /2 



4tt 2 B 2 A 



V8^a 2 f- 



2VttA 

3/2 



(C9) 
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where the definition of Q from Sec. HVBl was used. The 
correlation functions with dominant positional disorder 
are plotted in Fig. |U 



As discussed in the main body of the paper, the validity 
of this perturbative analysis is limited to short scales of 
the Larkin regime. To understand the behavior at longer 
distance, where the Larkin treatment breaks down, we 



utilize the renormalization group analysis. 

Appendix D: FRG analysis to the second order 

In this appendix we present the details of the FRG 
analysis by evaluating each term in Eq. (|83[) . Averaging 
over the high-wave- vector fields in the surface orienta- 
tional disorder term in Eq. (|83[) we find 




in which the cross term is linear in u > and thus vanishes 
in momentum-shell RG. The key part of the first term in 
the above equation is 

ii K [u% > (x,z)-u$ > (x,z)]^ 

Zn . 



"K fcfl 



(D2) 



in which Zq is the partition function of the high-wave 
vector components of the system 



Z>= [du^]e- H ^ T , 



(D3) 



and f al 3 is given by 

faP = C^ aa (0)-G> atP {0) 

f A d d - 2 q x f°° dq^T_ 
= 7/(1 - S aP )St. 



(i - <M 



(D4) 



For d — 3 it is easy to show that ?/ = 2 Jbx > wn i cn is 
a constant under the renormalization flow. We can see 
that with the expansion e~ K ^ al3 = 1 — k 2 f a p H , the 



zeroth order term contributes through rescaling and the 
first order term leads to the random oricntational pinning 
nonlincarity 



6Af 1 (u)=r ] A'}(u) 



(D5) 



with ' indicating derivation respect to u. 

Clearly, the average in the last term of Eq. (|D 1[) 
contributes to the surface positional pinning nonlincarity 
R v (u). To lowest order, we calculate the average 



l iK [ u S> ( x < z )- u o>( x ' z )} 



d, 



u 0>( x i z ) u 0>( x > z ) 



u o>( x , z ) - u 0>( x i z ) 



A d d - 2 q x 



l Ae -st (2^-2 

2C(1 ~ SaflW, 



oo 27T r q 



(1 - S aP ) 



(D6) 



while for d = 3 we have £ = 27T b\ > w hich contributes 

SRi la) (u) = -(A f (u)5l (D7) 

We now consider the surface positional disorder in 
Eq. (Sty, 
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a, 



^2T I ^ 2xdz { R v\ u o(x,z)-u^(x,z) 



a, p 



5^2T / 2xdz \ Rv u o<( x ' z ) ~ u oA x ' z ) 



^2T I dd 2xdz \ Rv u o<( x > z ) - u o<( x > z ) 



2 K 



2 R " 



u o<( x , z ) - u o<( x > z ) u 0> { x : z ) ~ u 0> ( x > z ) 
ito<(x,z) -u P 0< (x,z) ( u^ > (x,z)-w^ > (a;,z) 



(D8) 



in which the first term gives a correction through 
rescaling, the linear order term in makes a vanishing 



contribution, and the last term 

=»;<(u)M, (D9) 

[together with Eq. (|D7[l ] gives a first order contribution 
to the surface positional disorder. 

To study the second order contribution in the RG, we 
start with the (R v ) 2 term in (|55|) and obtain 




in which the average is given by 



u o>( x , z ) ~ u o>( x , z ) u o>( x ', z ') ~ Uol(x',z') 



2 r 



= [C%(6x, Sz)] 2 (S aia2 +6 aia2 6 l3l p 2 

-<W 2 <W 2 ~ <5a 2 /3i<W 2 )j ( D11 ) 

where Sx — x — x' and Sz = z — z' , and the three replica 
correction contributed by S aia . 2 term is irrelevant relative 
to the two-replica terms and is thus neglected. Then 
the second order contribution from the random positional 
pinning becomes 
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2T 



[C$(5x,6z)] \R" u^ k (x, z) - u^x, z) R" Uo<(a;', z') - Uq k {x' , z') 



-2i? 



— F 



£ <(0) 

[C^(fe,(5z)] 2 <^ R" Uq^x, z) - u% K (x, z) R" Uq k (x, z) - u^ K {x, z) 



+ 2 R " [ u o<( x ' z ) ~ u o<( x > z ) K" u o<( x , z ) ~ u o<( x > z ) 9 x Uq k {x,z) -u^ k (x,z) 



-2R'„: 



u^ K (x, z) - u^ K (x, z) i?"(0) 



(D12) 



With 



AS£ =^2 [C${Sx,Sz)]' 

•J 5x ,5 z 



d° 1 z q x f°° dq z T 2 



2tt r 2 



and 



Ao6£ 



2T 2 ./ Ae -« (2tt) 

(^) 2 [C>(fe,5z)] 



(D13) 



I 

pinning as 

6R<f a) (u) =A 
and 

SAf a) (u) = -A 2 R^(u)R^\u)S£. 



1 



4T 2 J5: 
4T 2 



1 f A d d -\ x r dq z T 2 q~ 2 



For c? = 3, we have 



7T - 2 



(D15) 
(D16) 

(D17) 



2n T 2 



(D14) 



8tt 2 J3 2 A 3 A 3 ' 



with only the main part of contribution, then we have The second order contribution of the random orienta- 
the second order contribution from the random positional tional pinning Af in (|83|) is 

I 



2t(^/: 



1 



32T 3 



E 



:c,z J x'z' 



Wo^x, z) - Mg^x, z) ^ Uq 1 (x, z) — Uq 1 (x, z) 



u^{x',z')-u^(x' lZ ') d x u^(x',z')-u p 2 (x',z') 



r 



(D19) 



which after separating the high and low wavevector 
components becomes 



d x 'U 



a-202 



A f (u^)\d x u^\ A, 

= A/ A, ((d x u^Y (d x ,u^ 



+2A f 



A'; <f 



+4A', A} («S^ 2 ) dtu^d^vg? 3 (uZ Pl uZ h d x u^d x ,u«f< 



(D20) 



where for simplicity of notation we defined [ug^x, z) — is given by (neglecting the three replica contribution) 
Uq 1 (x, z)] — 1/31 . The average in the first term of (|D20[) 
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The average in the second term of (ID20I) is 



[d x u^) 2 (d x ,u%f S 



= 4 



C^ g (Sx,5z) {8 ai a 2 5fi 1 2 - 2S ai p 2 S l 3 ll 3 2 ) . 



(D21) 



in which the correlation of the x derivative of the layer 



fluctuation, a , is 



d (8x, Sz) = (d x uo> (x, z)d x u (x + 5x, z + Sz)) > . 

(D22) 

Similar to (|D12|) . this term provides second order 
contribution as 



-A f (u)A f (u) - A/(«)A/(0) 



6RW(u) = A 3 
6Af b) (u) = ~A 4 A f (u)A'}(u)S£, 



se, 

(D23) 
(D24) 



where 



A.se 



2T 2 /, 



8x,5z 



q{6x, Sz) 



(D25) 



and 



d (5x, 5z) 



(D26) 



-da^Sfafa), (D27) 
where the correlation is defined through 

Cy 3j , 2 (fe, Sz) = (uq > (x, z)d x u(x + Sx, z + Sz)) > . 

(D28) 

It is easy to see that this term contributes to the random 
orientational pinning nonlinearity 



SAf c) 



-At 



A f {u)A" f {u) - A f {u)A" f {Q) 



-A/(0)A'» 



with 



A 5 



C> g/2 (Sx,5z) 



J§x,Sz 

The average in the third term of (|D20|) is 

\u 0> u 0> o x u 0> o x 'U Q> 1^ 
= 4 I C^(Sx,Sz)C^ d (Sx,Sz) + C^ g/2 (Sx,Sz) 



(D29) 



(D30) 



x{8 ai0l2 8p 1 p 2 —2d ai p 2 5p 1 p 2 ), (D31) 
which contributes to the random orientational pinning 
nonlinearity through 



SAf d \u) = -A e 



-A> f (u)A> f (u) - A> f (u)A' f (0) 



with 



(D32) 



A 6 - ft 



Sx.Sz 



C^(6x,Sz)C^ d (Sx,Sz) + g , 2 (5x,5z) 

(D33) 

The contribution to the Hamiltonian from the cross 
term, the third term in (ID20[) , is 



1 (Ha.H^I = J2 



T 



x.z J x'z 



u^{x,z) - u^(x,z) d x Uq 1 {x,z) - u^(x,z) 



xR, 



[u?(x',z')-ufr(x',z') 



(D34) 



Similarly, we can separate the long and short scale field 
and keep only the leading order contribution, thus 
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r 



The average in the first term of (|D35[) is easy to calculate 

= 4 [C£ (Sx, Sz)] 2 (S aia2 S Pl p 2 - S ai/ 3 2 S0 1 p 2 

and it contributes 

6Af e \u) = A\A'}(u)R':(u) ~ A';(0)<M 



-A'»<(0) 



51, 



(D37) 



with constant A given in (|D13[) . 

The average in the second term of (|D35[) is 



_ ^2ft^ife), (D38) 
and it contributes 



= -A, 



A f (u)R':(u) - A f (0)I%(u) 



-A/(«)<(0) 



S£, 



(D39) 



with Ac, given in (|D30p . 

These second and lower order contributions sum 
up into flow equations of the non-linearities given as 
Eqs. flU and (JHU) in Sec. IVAl 



Appendix E: Full behavior of long scale correlation 
function 

The RG and matching analysis in Sec. [V] predicts that 
the long scale smectic phonon correlation is dominated 
by random orientational pinning, with strength A j (T) = 
AfSf(T). The correlation is given by 



C(x,z) = 2 



dq z dq x A f (T)q 2 x \J\ 2 q x + q 2 - Xq 2 x 



27T 27T 2B 2 X q 2 (X 2 qi + q 2 ) 

x (1 - e iq * x+iq * z ) , (El) 
in which the integrations can be safely extended to 
infinity with our focus on long scale behavior. Given 



A as a positive number, it is easy to observe a relation 
that 



C{VAx,Az) = VAC(x,z), 



(E2) 



which is consistent with the results shown in Eq. (JT)), on 
the x or z axis 

C(x, 0) oc x, and, C(0, z) cx ^fz. (E3) 

The correlation can thus be characterized by 

' \z s 



C(x, z) oc xf 



in which at different limits of the ratio t = Xz/x 2 



f(t) = const., for t — >■ 0, 
f(t) oc \/t, for t — > co. 



(E4) 



(E5) 



Redefining variables 

(u,v) = (x 2 ql,zq z ), (E6) 
so that (q x ,Qz) = {\/u/x,v/z), Eq. (|E1[) then becomes 



C ^ = ^xJ d «*j ^ q 2 (X 2 q x + q 2 ) 
x [1 -cos(q x x) cos(q z z)] 

A ^ x 2 t 2 r du r w^^-^ 



4tt 2 b 2 a 3 "^ y -■"./„ ~" v ~ v 2 (t 2 u 2 + v 2 ) 

x [l — cos (Vu) cos (v)] 
A/(T) 



4tt 2 B 2 A 3 



s/(f), 



(E7) 



with 



oo poo 



f(t) = 2r / du / dwV" 



-\ft 2 U 2 + f 2 — tit 
1,2 f +2 7/ 2 _i_ t ,2 ^ 



IQ JO V 2 (t 2 U 2 + V 2 ) 

x [l - cos (y/u) cos (w)] . (E8) 

Writing 1— cos (\/u) cos (v) as [1 — cos (y / u)]+cos (y/u){l — 
cos(i>)], it is easy to obtain the full analytic expression 
and limits of f(t) as given in Eq. (|104[) . and plotted in 

Fig.ra 
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